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Asymmetric current -voltage {I{V)) curves, known as the diode or rectification eftect, in one- 
dimensional electronic conductors can have their origin from scattering off a single asymmetric 
impurity in the system. We investigate this effect in the framework of the Tomonaga-Luttinger 
model for electrons with spin. We show that electron interactions strongly enhance the diode effect 
and lead to a pronounced current rectification even if the impurity potential is weak. For strongly 
interacting electrons and not too small voltages, the rectification current, Ir = [liV) + ^(~V')]/2, 
measuring the asymmetry in the current-voltage curve, has a power-law dependence on the voltage 
with a negative exponent, Ir ~ V~^^\ leading to a bump in the current-voltage curve. 
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I. INTRODUCTION 

Current rectification, also known as the diode effect, 
occurs whenever the transport characteristics in elec- 
tronic conductors are asymmetric under the application 
of the driving voltage. This asymmetry can have various 
origins. Best known perhaps is the mismatch of band 
structures in macroscopic diodes, which leads to a po- 
tential wall that blocks the motion of the particles in one 
direction, while it is not seen by the particles moving in 
the opposite direction. 

In the last years, rectification in mesoscopic and single- 
molecule systems has attracted much attention. Even 
though rectification by single asymmetric molecules has 
been suggested 30 years agoi, it has been realized ex- 
perimentally in the 1990s onljiS*^. Asymmetric wave 
guides were constructed in the inversion layer of semi- 
conductor heterostructuresi^. Transport asymmetries 
have been observed in carbon nanotubesiS and in tun- 
neling in quantum Hall edge states^,. The experimen- 
tal progress has been accompanied by much theoretical 
activitjiifliiLiSiiiiiiiSiiSiil with the main focus on Fermi- 
liquid systems. 

An interesting source for rectification arises when 
the particle interaction is strong enough that a single- 
particle description becomes invalid. This naturally oc- 
curs in one-dimensional systems, for which a Luttinger 
liquid behavior is expected. It has been known for some 
time now^ft that in such systems the current is strongly 
affected by the presence of impurities, even when they 
are weak, due to their renormalization by the electron- 
electron interactions. One can thus expect that asym- 
metries in the impurity distribution leads to strongly 
asymmetric current-voltage curves. This issue was ad- 
dressed very recently in the framework of Luttinger liq- 
uids of spinless particlesii. The rectification current 
Ir{V) = [I{V) + I{-Vy\/2 can be measured as the dc 
response to a low-frequency square voltage wave of am- 
plitude V and expresses the asymmetry of the I—V curve 
for forward and reverse bias. It was shown that a single 
weak asymmetric impurity is sufficient for a pronounced 
rectification effect (see Fig. for a sketch of the system), 
leading to a large rectification current, /r(^), at low volt- 



ages V . Moreover an unusual behavior of the current was 
revealed for systems with strong repulsive interactions: A 
power-law dependence of /,■ on the voltage with a nega- 
tive exponent, Ir ~ T^"'^', within a range V* <V < V** 
[with V* and V** being expressed by some powers of 
the impurity potential U ; see also below] , i.e. the recti- 
fication current increases as the voltage is lowered. At 
V <V* the increase crosses over into a regular decrease 
such that the equilibrium condition / = at = is 
met. The qualitative behavior is shown in Fig. |2| 



In this paper we pursue this work for real electrons 
carrying a spin. We show that we can provide lower and 
upper limits for the voltage V , within which perturbation 
theory can be applied and determine the leading power- 
law behavior of the rectification current, Ir ~ , as 
a function of the charge and spin interaction strengths, 
Qc and gs- This leads to a phase diagram for the leading 
power-law dependence of Ir with a much richer structure 
than in the spinless case (Fig. 0). We show that, similar 
to the spinless case, a negative exponent z, Ir ~ F"'^', 
appears for strong repulsive electron-electron interac- 
tions within a voltage range that is determined by the 
bare impurity scattering strength. Fig.|3shows our main 
result, the leading power-law behavior of the rectifica- 
tion current Ir- The hatched region in the figure marks 
the range of (gc,gs), in which the exponent z becomes 
negative and in which the rectification current shows the 
behavior sketched in Fig. |21 



The paper is organized as follows: In the next section, 
we introduce the model of electrons in a one-dimensional 
system, and argue how scattering on a single impurity 
can lead to rectification. We then quantitatively address 
this problem within the bosonization approach, and show 
that the most relevant power-law behavior of the rectifi- 
cation current can be obtained within second and third 
order perturbation theory. The results are listed in Ta- 
bleland in Fig. 121 In the appendix we show that higher 
order perturbative contributions cannot exceed the lead- 
ing second and third order expressions. 
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II. MODEL AND ORIGIN OF RECTIFICATION 

Let us consider a one-dimensional conductor of length 
L with electron-electron interactions that are effectively 
short-ranged due to screening by gates. Such a system 
allows a description by the Tomonaga-Luttinger model, 
given by the Hamiltonian 



H 



Ax 



n- 



-hvi 



dxdy^K„^>{x - y)il;l{x)il;l,{y)-ipa'iy)i^a{x), 

(1) 



where ipRa and t/'Lo- are the operators for right- and left- 
moving electrons with spin a |, and ipa = ipua +4'Lct 
is the conventional electron operator. Uaix) is the asym- 
metric potential (i.e. U{x) ^ U{—x)) localized in a small 
region about x = 0. K{x — y) describes the electron- 
electron interaction. For the following discussion it is im- 
portant to assume that the long-ranged Coulomb inter- 
action is screened by the gates, so that K{x — y) becomes 
a short -ranged, rapidly decreasing function of {x — y). 

On its two ends, at a; = ±L/2, the system is adiabati- 
cally coupled to electrodes that serve as reservoirs of par- 
ticles, and whose chemical potentials /ii^2 are controlled 
experimentally: We assume that one electrode is fixed at 
the ground, /i2 = /i, while the other one is connected to 
the voltage source, fii = fi -\- eV. For such situations, it 
is possible to distinguish between two effects, addressed 
more quantitatively below, that lead to rectificationii: 

(1) The "injected- density driven" rectification as the 
result of the dependence of the charge density on the volt- 
age: For simplicity, let us consider noninteracting parti- 
cles. For the voltages ±V, only electrons in the energy 
ranges between [/i, fi-\-eV] and [fi—eV, ^] (at zero temper- 
ature) can contribute to the transport. The presence of a 
scatterer U with an energy-dependent transmission coef- 
ficient R{E) in the system leads to different transmission 
coefficients for ±V and thus to different currents. This is 
seen as follows: For noninteracting particles, the reflec- 
tion coefficient R{E) is independent of the propagation 
directionia, and the current is I{V) - [1 - R{E)]AE. 
If the bandwidth Ep is the only relevant scale for the 
energy dependence of R{E), then, for small U and V, we 
have R{E) ^ U'^/Ep and the rectification current Ir ~ 

J^^[R{EF~E)~R{EF + E)]dE w 2R'{Ef) Jq^ EdE - 
R(EF){eV)'^/EF ~ U^{eV)yEj,. The rectification cur- 
rent is nonzero. For noninteracting particles it is pro- 
portional to V^. As shown below, a modified power-law 
dependence on V is obtained in the presence of electron- 
electron interactions. 

We remark that this argument does not require a spa- 
tially asymmetric impurity potential since the asymme- 
try is introduced through the injected charge densities. 



(2) The "asymmetric-impurity driven" rectification ef- 
fect is independent of injected densities. It results from 
the renormalization of the asymmetric potential U{x) 
by the electron-electron interactions, which leads to the 
asymmetric current-voltage curves. For the Luttinger 
liquid, this naturally involves multi-particle processes, 
so that all such possible terms have to be taken into ac- 
count in the modeling. As shown in Ref. 0, this effect is 
absent in the first two orders in the scattering potential 
U, and we must address it perturbatively at the order 

The "asymmetry-driven" rectification effect can be 
qualitatively visualized (in a mean-field picture) as fol- 
lows: In an interacting system, electrons are backscat- 
tered by a combined potential U{x) — U{x) -\- W{x), 
where W{x) is the self-consistent electrostatic potential 
created by the average local (nonuniform) charge density. 
Depending on the voltage sign, the density decreases or 
increases as a function of the position x and the magni- 
tude of U{x) (see Fig.P). Asymmetric impurities create 
different density profiles for opposite voltages, and lead to 
different combined backscattering potentials U{x). The 
rectification effect is a consequence of the modification of 
the current by the backscattering potential. 



III. BOSONIZATION 

For the quantitative treatment of the rectification ef- 
fect we use the bosonization technique, which has be- 
come a standard tool for one-dimensional problems^Si. 
In the bosonization language, the creation of a right- or 
left-moving electron at the coordinate x is expressed by 
ipt^{x) ^ r]l^e-Kp{-\-i[±kF^x -\- (p^aix)]), for ly = R,L, 
where kFu/T^ are the densities of right- and left-moving 
particles, where 77^^, the Klein factors, raise the total 
number of right- or left-moving electrons with spin a 
by one, and where ip^aix) is a boson field that describes 
the dressing of this particle by a chain of particle-hole 
excitations. In correlation functions, the Klein factors 
keep track of particle conservation and can contribute 
with signs to the expressions, as they anticommute for 
different {v^a). These signs are all equal in the present 
calculations since we assume conservation of the Sz com- 
ponent of the spin by the scattering process. This allows 
us to drop the Klein factors in the following expressions. 

In one dimension, the charge and spin degrees of free- 
dom of the Hamiltonian decouple, and it is con- 
venient to set (pc — J2u=L fll'^i'T + and (ps — 
J2iy=L rI'I'i'^ ~ 'Pi^l]/'^! which are bosonic fields related 
to the charge and spin densities as pc = e[dx<t>c + 2fci?]/7r 
and ps — {h/2)dx<t's/T^, where kp — [kpR -\- kp]f\/2. 

For U = Q the bosonization leads to a quadratic action 
for the fields (j)c s , which can be written in the form 



so=j dxdt 8^ m^a) 



(2) 
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The quantities and Qs are the charge and spin in- 
teraction strengths resulting from the screened electron- 
electron interaction K{x — y). A noninteracting system 
is characterized by = 5s = 1/2. Repulsive (attrac- 
tive) interactions are expressed hy gc < 1/2 {gc > 1/2). 
Interactions with conserved SU{2) spin symmetry have 
gs = 1/2. For gs < 1/2, a neglected sine-Gordon term 
of 0s in Sq, describing spin-flip backscattering, would 
become relevant. The physics of such systems are con- 
siderably different to the conductors described here (see 
e.g. Ref. l20|) . since the fields (j)s would order and become 
massive. We exclude such situations explicitly by assum- 
ing gs > 1/2, such that all backscattering terms in iSo arc 
scaled to zero. 

If we assume that U{x) is concentrated in a small 
region ajj <^ L about a; = 0, the backscattering can 
be described by the values of the fields at a; = only, 
$c,s(i) = (f'c.six — 0,t). This "single-point" description 
of U{x), restricts the model to wavelengths longer than 
au, i.e. to energy scales smaller than Eu <^ hvp/aij 
(with vp the Fermi velocity). Since the characteristic 
energy scale is set by the voltage, Ejj provides an upper 
limit for the applied voltage, eV <C Ejj, and the validity 
of the model. Under this assumption, the impurity term 
has the formic. 



Si^p = ^J dt C/(ne,n,)e*"("-"=)e™=*=+™=*^ 

7ic ,ns 

(3) 

The integers ric^s physically express the transfer of Uc 
charges and rig spins between the right- and the left- 
moving electrons (where, for instance, ric — Ug = 1 cor- 
responds to the backscattering of an up-spin electron, 
and Uc = 0, rig = 2 to the backscattering of an up-spin 
electron and a down-spin electron with opposite inci- 
dent directions). U{nc,ns) and a{nc^ns) are modulus 
and phase of the effective multi-particle backscattering 
potential (for symmetric potentials, U{x) = U{—x), the 
a(nc,rts) would vanish). Forward scattering can be ab- 
sorbed in Sq by a shift of (pc- Since charges and spins 
are bound to physical electrons, the sum (indicated by 
the ' next to it) runs over nc,s such that ric + is even. 
Iricl = I'^sl = 1 corresponds to the usual 2kp backscat- 
tering of electrons. The most relevant contribution to 
the currents arises from the backscattering processes with 
ric + ns = 0, ±2. These terms dominate in the "density- 
driven" rectification effect, but higher orders in ric + i^s 
are required for the "asymmetry-driven" rectification ef- 
fect, and the full sum over n^s, therefore, must be kept. 

Throughout the following analysis, U{x) is assumed to 
be weak, |?7| ^ Ep (with Ep ~ bandwidth). Due to the 
strong renormalization of the impurity strength through 
the electron-electron interactions, the precise nature of 
U is of no major importance; U{nc,ns) is an effective 
many-particle potential, dressed by short-time electron 
processes with frequencies well above Eij. Its magnitude 



can be roughly estimated asi^ 

U{nc,ns) r^kp j dxe2''=^"'^^C/(x), (4) 

which means that U{nc,ns) ~ U. Its further renormal- 
ization by low-energy processes, leading to a power-law 
correction in V, is considered below. The coefficients 
U {ric, Us) depend further on the applied voltage through 
the densities kpL and kpn. 

Since the Hamiltonian is Hermitian, U{nc,ns) — 
U{—nc,—ns) and a{nc,ns) — ~a{~nc,~ns). Because 
of the spin symmetry, U{nc,ns) ~ [/(ric,— rig) and 
a{nc,ns) = a(nc,—ns). For spatially asymmetric po- 
tentials, as considered here, a{nc,ns) ^ a{—nc,ris). 

If [/ = 0, the system is characterized by right -moving 
particles injected from the left reservoir with the chemi- 
cal potential /ii, and left-moving particles injected from 
the right reservoir with chemical potential fi2 ■ Due to the 
absence of backscattering in Eq. jSJl, the right- and left- 
movers are in equilibrium with the corresponding reser- 
voirs. In the presence of an impurity U, the equilibrium 
notion of chemical potentials becomes meaningless in the 
system; the quantities fj,i^2 enter only through the volt- 
age dependence of the injected carrier densities in the 
description and appear in the backscattering term that 
couples right- and left-moving particles. Explicitly, this 
can be seen by switching to the interaction representa- 
tion H ^ H = H - fiiNR - /iaiVi for U = 0. This 
shifts the single particle energies of the injected particles, 
and the electron operators become ipRi^i t) ~^ ''Pr{x, t) = 
e^'^i*/ Vfl,(a;, t) and ?Al(x, t) '4}l{x, t) = e*'"=*/VL(a;, t). 
Hence, for ?7 7^ 0, the backscattering of every charge de- 
scribed by (Simp acquires a time-dependence on the dif- 
ference of chemical potentials, /ii — /i2 — gV , as 3='=*^^*/''. 
For the ric backscattered charges, this becomes 

Sin.p^- fdt [/(n„n,)e^"("='"') 

•' na,n,eZ (5) 

The time-dependence of the backscattering operator re- 
flects the nonequilibrium constraints. In order to cal- 
culate expectation values, we have to consider an ap- 
propriate nonequilibrium technique, such as provided by 
the Keldysh methodi2i42^: We assume that in the re- 
mote past, t —00, the system is fully described by 
Sq and the ground state |0) of the shifted Hamiltonian 
H — fiiNfi — IJL2NL. Averages are taken over this well- 
defined ground state |0) only, while the impurity term 
5imp is taken into account through the time evolution 
operator S'(— oo,t), given by 

S{t', t") = Texp 1^-^ dtiJi„,p(t)^ , (6) 

with -ffimp the impurity Hamiltonian in the interaction 
representation, and T the time order operator. The ex- 
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pectation value of an operator 0{t) is expressed by 

{0{t)) = {0\S{~oo,t)O{t)S{t,-oo)\0), (7) 
with S{t, —oo) — S{—oo, ty . 



IV. CURRENTS 

In a clean system, the current flow is proportional to 
the voltage, and is given by the Landaucr formula^^'^^'^^, 
Io{V) = 2Ve'^/h (the factor 2 accounts for the two spin 
channels). In the presence of the impurity scatterer 
U, the backscattering corrects the current as I{V) = 
loiy) + Ibs{V)- This shows that the rectification current 
depends only on the backscattering currents, IriV) — 
[Ihsiy) + Ibs{—V)\/2. The backscattering current opera- 
tor can be obtained, for instance, by the time variation of 
the number Nji of right-moving particles, = ^Nji = 
i[H,NB\/h. If we set e = ft = t;F = 1 for the ease of 
notation, this yields, in the interaction representation. 



hs{V, t)^iJ2' nMinc, n,)e^"("-"=)e™=*=+™=*=e^ 



iricVt 



(8) 



and we need to calculate 



{hsiy)) = (0|5(-(X3, 0)/b,(F, 0)5(0, -oo)|0), (9) 

where S{t^ t') is the time evolution operator, and |0) the 
ground state of the system described by Sq. 



V. RESULTS FROM WEAK COUPLING 
THEORY 

For a weak impurity potential, the currents can be cal- 
culated within perturbation theory. 

As shown by Kane and Fisher—, the backscatter- 
ing potential is renormalized by the electron-electron 
interactions and scales as a power-law with the char- 
acteristic energy (here set by V) as Uctf{nc,ns) ^ 
U{nc,ns){V/ Epy^'^^-^^'"^^^^ ■ Perturbation theory is ap- 
plicable when Ucs <^ Ep. 

For n^gc + n^ds > 1, this condition is always met for 
y, {/ ^ Ep, while for n^gc + n^gs < 1, wc must have 



V ^V* 



(10) 



for any admissible values of ncUs- We set V*^ ^^^ — 
whenever n^gc + "-sffs > 1- Since > 1/2, only Us — 
0, 1 have to be considered. We always have 



2n+2,0 



< 



^2*n,o ^"^d V^n+i^i < ^2n-i,i5 that the only important 
quantities are V-^i and V2q- For gc + gs < 1 and < 1/4 
both and V^q are nonzero, and we have V{i > V^q 
for gs < 3gc- The lower cutoff energy is given by 



V* = max[yi*i,y2*o] 



(11) 



Interestingly, if V* > 0, perturbation theory requires a 
voltage V that is not too small, V > V* , which is much 
in contrast to the usual perturbation theory, in which V 
would be required to be a small parameter close to = 0. 

The upper energy limit is set by the energy Eu, at 
which the corrections to the model become important. 
The perturbation theory is valid in the range eV* <C 
eV^Eu- 

In this case, the rectification current is dominated by 
the second and third order perturbative expressions from 
an expansion in powers of U. Due to the renormalization 
of the potentials, and due to the constraints on particle 
conservation and even sums of ric + rig, third order ex- 
pressions can become larger than the second order ones 
and have to be taken into account. In the appendix we 
further show that, for V > V*^ higher orders cannot ex- 
ceed the second and third order expressions, and can be 
safely neglected. 



A. Second order 

As mentioned, the rectification effect arises from the 
asymmetry of the charge/spin density profiles. The 
"injected-density-driven" rectification effect is due to the 
dependence of the density on the voltage, and hence de- 
pends on how the coefficients U{nc,ns) change upon the 
variation of the density (see Eqs. Q and l@}). The dom- 
inant contribution appears at second order in U, which 
reads after expanding Eq. @ 

where we have used the notations n — (nc,ns) and 
n^{t) — fic^dt) + Us^sit)- (•■•)o is the average over 
the ground state of Sq, Ck is the Keldysh contour 
— (X) — > — > — cx), and Tc is the chronological order op- 
erator on Ck- Particle conservation imposes m = ~n. 
Since a{—n) = —a{n) the phases a cancel each other. 
Ignoring the voltage dependence of U, this expression, 
therefore, changes the sign upon V ^ —V, and would 
naively not contribute to the rectification effect (which 
is a consequence of the invariance of the action Sq under 
the change (j)c.s —4>c,s)- However, since U depends on 
the voltage through the density kp (see Eq. (@J)), we can 
expand it in powers of F as t/ = t/o + VUi -I- . . . . At 
small voltages, V <C Ep, the correction is linear, and if 
the bandwidth Ep is the only relevant scale for the en- 
ergy dependence of C/, we have Ui ~ ^ Ua/Ep. The 
rectification current can be written as (choosing > 0) 

ll^\V) = 2V I dt^'in,t/o(n)C/i(n)e-™=^* 

X ^Tce™*^'^) e-^"*(*)^^. (13) 
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The propagators in the latter equation have the usual 
form known from the Luttinger liquid theorj*2£: 

^rp gine*c(0) + ins*s(0) g-mc*c(t)-™s *s (*) ^ 

= gn^(T,*,(0)<I.e(t))o+n'{Tc*s(0)<E's(t)>o^ (14) 

where the free propagators are given by 

($c,s(ti)$c,s(i2))o = -2.gc,.ln(i(ii - h) + 6), (15) 

with (5 > an infinitesimal constant. We can now extract 
the voltage dependence of the current by changing to the 
time variable s ~ Vt. The propagators become 

/ rp mc*c(0)+m,*,(0) -inc*c(t)-ins*s(t)\ 

^Jee e 

From the integration measure, we obtain another factor 
F^^, which cancels the V from the expansion of the po- 
tential. Hence the rectification current has the form 

^^'H^) = I]'"cC(n)y""'(") (17) 

with C constants of the order of C ^ ?7oC^i • Ij and 

z(2)(n) = z(2)(nc,n,) = 2nlg, -f 2nlgs. (18) 

The leading order is obtained by minimizing the expo- 
nent, which is achieved by one of the processes (nc, 't-s) = 
(2, 0), (0, 2), (1, 1) (and the combinations obtained by 
changing signs) only, as backscattering of more particles 
leads to larger z'^^^. Terms with Uc = do not cou- 
ple to the voltage (e™<=^* = 1) and coincide with the 
equilibrium value at 1^ = 0. Therefore, their amplitude 
vanishes. The comparison of the remaining two processes 
shows the dominance of {nc,ns) = (1, 1) [labeled as (A) 
in Fig.E] for gs < 5gc and (nc,ns) = (2,0) for gs > 3gc 
[(B) in Fig. 121. The noninteracting system is character- 
ized by = 5s = 1/2, leading to z = 2. This result 
certainly agrees with the result from the application of 
the Landauer formalism to a noninteracting problem. 

B. Third order 

The "asymmetry-driven" rectification effect appears at 
third order. The contribution to the backscattering cur- 
rents becomes 

X t/(n)t/(TO)C/(/)e^"(")+'"(™)+'"(') 

X ^Tce"*(°) e""*(*i) e''*(*^)^^, (19) 

which has to be evaluated with the constraints -I- -I- 
Ic = Us + rris + Is = 0, together with ric + Ug, rric + ttLs, 



and Ic + Is being even. Again, diagrams in which ric = 
~ Ic ~ 0, do not couple to the voltage and vanish. 
The third order term is no longer invariant under the 
transformation {V, $c) ^ (— V^, — <i>c) because, generally, 
a{n) + a{m) + a{l) ^ (i.e. a{nc,ns) + a{mc,ms) ^ 
a{nc + mc,ns + rus))', the invariance would require the 
additional transformation a — *■ —a, corresponding to a 
mirror refiection x — > —x of U . Since the sum of the ex- 
ponents a does not vanish, /,. ^ 0, and the same analysis 
as before leads to 

7^3) (y ) ^ c{n, m, l)V'"^ , (20) 

with 

z^^\n,m,l) = g,{nl+ml+ll)+gs{nl+ml+f,)-2, (21) 

and the amplitudes C ^ U^. Tabled lists the smallest 
exponents z at this order, characterized by the letters (C) 
and (D), for which the rectification current is nonzero. 
The exponent zo is smaller than zc for gs > 9gc- 

To obtain the leading power-law behavior for given 
{gc,gs), we compare the amplitudes from second order, 
^ and third order, U'^V^'' \ for voltages 

V ^ V* using Eq. H10|) . If V* = 0, the mere compar- 
ison of exponents is sufficient. The result of comparison 
is represented in Fig. O The diagram is not a phase dia- 
gram in the proper sense because the crossover lines shift 
with the voltage. Close to the crossover lines, all con- 
tributions of the neighboring phases are large, whereas 
far from the boundaries, a single power-law dominates. 
Finally, as shown in the appendix, higher order pertur- 
bative corrections lead to less relevant power-laws and 
can be neglected. 

VI. CONCLUSIONS 

The above results show that the inclusion of spin de- 
grees of freedom leads to a more diverse behavior of 
the rectification current for different interaction strengths 
than in the spinless case of Ref. The leading power- 
laws for the rectification current as a function of the in- 
teraction strengths gc and gs are represented in Fig. |31 

The most interesting region is characterized by a nega- 
tive exponent, z < 0. It leads to the unusual behavior of 
a decreasing rectification current as the voltage V > V* 
is raised. Since negative exponents appear only from the 
third order contributions, they are a realization of the 
"asymmetry-driven" rectification effect, and due to the 
presence of an asymmetric impurity in the system. The 
decrease stops at the upper voltage V** ~ {7i/(l^*^'l+^*^')^ 
where the amplitude of the second order contribution ex- 
ceeds that of the third order. The qualitative behavior is 
sketched in Fig. [3 

The rectification effect is strong if the magnitude of 
Ir becomes comparable to that of the total current 
I{V) ^ V and to that of the most relevant contribution 
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to the backscattering currenliia, Ibs{V) ~ y2n^3,+2n^g,-i 
with (jicrig) = (1,1) or (nc,ns) = (2,0). For interaction 
strengths gc,s such that V* = 0, the corrections are weak 
and generally IriV) < Ibs{V) < I{V). For V* > 0, 
however, Ibs{V) becomes comparable to I{V) at V ^ V* , 
i.e. close to the limit of validity of perturbation theory. 
For such voltages, in regions (C) and (D) the ratio of 
rectification current to total current can roughly be es- 
timated as Ir{V)/I{V) - [/3(y*)^*"-i. In phase (C), 
when V* — V{i (i.e. 3gc — 5s > 0), this leads to the 
ratio Ir{V)/I{V) - (V*f3c~gs ^ jj{3g.-a.)/(i-a.-9,) ^ 
The rectification current, therefore, can become compa- 
rable to I{V) (as well as to Ibs{V)) when ig^ — gs is 
small. The comparison in the region where V* = V^q (i.e. 
3(7c — <?s < 0) yields a similar condition for 2gs — 6gc be- 
ing smah. In phase (D), we find Ir{V)/I{V) - (V"2to)^^^% 
which may lead to an enhancement of the rectification ef- 
fect for weak impurities as gc becomes small. 

To conclude, we have shown that one-dimensional 
electronic conductors with screened electron-electron in- 
teractions can exhibit a strong rectification effect in 
the presence of spatially asymmetric impurity scatter- 
ers. The rectification current shows a power-law behav- 
ior, Ir ~ V^, with an exponent z that depends on the 
electron interaction strengths gc and gs only, and which 
can be determined from second and third order perturba- 
tion theory. For small values of gc the exponent becomes 
negative, leading to the upturn of Ir about the voltage 
V* (Fig. I2J). We note that this effect was obtained in 
the framework of weak coupling theory. For voltages 

V < V* , the effective impurity scattering strength ex- 
ceeds Ep, marking a crossover into the strong backscat- 
tering regime^^'^^. All currents eventually vanish at 

V = 0. 
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APPENDIX A: ESTIMATE OF HIGHER 
PERTURBATIVE ORDERS 

In the preceding study we have shown that third order 
contributions can exceed the second order ones for small 
values of gc,gs- To complete this, we have to give an 
additional proof that higher perturbative orders iV > 4 
cannot exceed these values in the physical range of > 
and gs > 1/2. 

Since the calculation of the rectification current in- 
volves contributions that are lesser relevant than the 
leading power-laws for the backscattering current, we 
have also to show that the neglected backscattering 
term in So, proportional to J dx cos{2(j)s), cannot gen- 
erate additional important corrections to the rectifica- 
tion current. Under renormalization, this term has 



to be completed by allowing the general interaction ~ 
J dx cos{2'hs(t>s) with a summation over integer hs 
(the notation fis is used to distinguish these indices from 
those Us appearing in U{nc, ««)). Under the renormaliza- 
tion group, when integrating over high energies down to 
the characteristic energy V , such terms are renormalized 
asSa (l^/£;^)(2"»)'9=-2. 

We choose the following strategy: Let us assume to be 
in the perturbative region with voltages V > V*, such 
that U,V <. Ef and C/(V"/£;f)»="'+»»"'-i < Ep. A 
general higher order correction to the current has the 
form 

{nc,ris} {"^s} 
A'^ factors M factors 

(Al) 

where for simplicity we set Ep — 1 and neglect pref ac- 
tors of order 1. Charge and spin conservation here re- 
quires that ^ ric = and X^^s + Sl^'^s) = 0. Since 
gs > 1/2, the product over the ~^ is always < 1, 

and if we denote by /(^) the expression obtained from 
Eq. ^ by suppressing these M factors, we always have 
j(N,M) <- j(N) ^ Therefore, the corrections arising from 
/dx cos (2 0s) are always small. 

Let us choose a subset of factors of and denote 

this quantity by J, 

j = v Yl (;7y'''="'+^'="=-i). (A2) 

selection of {nc,ns} 

Since all factors arc smaller than 1, we have 

j(NM) < j(N) ^ J (A3) 

We will now show that, for > 4, it is always possible 
to choose a J that is smaller or equal to the dominating 
second or third order contribution. The larger the nc,s, 
the smaller become the factors in /(^' and J. Much in 
the analysis, therefore, depends on the maximal value of 
I Tic I, which we denote by fie- 

This method of comparison does not hold for the com- 
parison between the second and third order expressions 
themselves. Since these expressions are products of two 
or three factors only, we cannot choose a suitable subset 
for J but have to deal with the full expressions. The 
exponents are therefore largely determined by the con- 
straints of particle conservation plus having even sums of 
charge and spin numbers. The result of comparison was 
discussed above. For A^ > 4, however, the larger freedom 
of choice of J allows us to show that these higher order 
expressions are small compared to the second and third 
order ones. The following estimates contain the proof of 
this statement. 

For any choice of maximal fie we choose a suitable J 
consisting of two or three factors [/T/"cffc+"sffs-i ^jjg 
prefactor V or, if required, the prefactor with the 
additional V arising from the expansion of a J7 (similarly 
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to the second order calculation) . Wc then show that this 
J is smaller or equal to one of the expressions 

= [/2y2f<=+2s=, (A4) 
Ib^U^V'^S", (A5) 

for any 5c > and Qs > 1/2. This will prove that 
/^^^ < J cannot exceed the dominant second or third 
order expression, given by the maximum of Ia,b.c,d- 

Since fic > 1, the particle conservation imposes that 
another or several numbers are different from zero. 
These are denoted by he and rhc below. If not further 
specified we only assume that \fic\, \ fhc\ > 1- If fic is odd 
there is at least another odd \t%c\ > 1 and, due to the 
requirement of even ric + n^, there exist at least two odd 
\ns\ > 1, denoted by hg and fhs- 

We distinguish between the following five cases: 

1) nc > 2 is even; all nonzero \nc\ = fic'. If in the couples 
(nc,ns) all Uc ^ Q have = 0, /(^) (and 

is symmetric under the change of sign V ^ —V, and 
the rectification current vanishes (note that the rig for 
which ric = and the fig can be nonzero though). This 
situation is similar to the second order case discussed 
above. A rectification current exists since the poten- 
tial U depends on the voltage through the density kp- 
An expansion of a to linear order in V allows us to 
choose a J of the form 

J = V'^{UV^cac-'^){UV^o9c-l) < [/2y(2=+2=)ge ^ 

(A8) 

On the other hand, if not all are zero, there must 
be an even = hs with \hs\ > 2, and we can choose 
for J 

J = y(t/y«c<'c+n?g.-l^^^yfi^9e-l) 

< [/2y(2^+2^)flc+2=sa-l < u2y8gc+2gs < j^^^^ (^9) 

where we have used gs > 1/2. 

2) fic > 2 is even; all nonzero \nc\ are even but are not 
all equal: This condition excludes fic = 2 as it coin- 
cides with the previous case. For fic > 4 there exist 



two other nonzero and even \nc\, \ific\ > 2 or another 
\hc\ > 4 and an arbitrary \rhc\ > 0. We then choose 



3) nc > 2 is even; there is an odd |nc|: Since X^ric = 
there exist at least two odd jnd, |toc| > 1, and 
since ric + rig must be even, there must be two odd 
I'^sli l^i'sl > 1- This leads to the bound 

J = [/30"c+"c+'"c)9c-(-(fte+™a)9s-2 

< [/3^(2^ + l + l)ffc + (l + l)3.-2 ^ ^A^^) 

4) fic = 1: In this case all ric are \nc\ = 1 or zero. If all 
nonzero ric have Us = ±1, l'^^^^^ is symmetric under 

y — > — and wc have to expand a [/ to linear order 
in V. If we choose a |fic| = 1 and two \hs\, \ms\ = 1, 
we can set 

J = y^2y(nl+nl)gc + {f>''l + ml)g,-l 

< f/2y(i+i)sc+(i+i)s. = (A12) 

On the other hand, if there is a couple {fichs) with 
fic 7^ and |fis| > 3, we have the estimate 

J = [/2y("c+»c)Sc-FftJs.-l < [/2y(l + l)3,+32g,-l 



5) fic > 3 is odd: Since fic is odd there must be at least 
another odd |fic| > 1, and, to fulfill that ric + ng is 
even, there must be two odd |ns|, |ms| > 1, such that 

< j72y(3^ + l)9c + (l + l)9,-l < [/2ylOSc < J^. (A14) 



We conclude that higher order terms cannot exceed the 
second and third order expressions for V ^V*. 
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FIG. 2: Qualitative behavior of the rectification current, 
Ir{V) = [I{V) + I{-V)]/2, in the region, where the lead- 
ing (third order) exponent is negative, z < 0. For voltages 
larger than V* , the rectification current decreases until, at 
V** , the leading second order term exceeds the third order 
amplitude, and leads to a further increase of the current. 



FIG. 1: Sketch of system of length L with a simple asymmet- 
ric potential U{x), consisting in two point scatterers of differ- 
ent magnitude in a region au <S L. The horizontal lines qual- 
itatively represent the density profile (averaged over Friedel 
oscillations) expected in the system. The reservoirs determine 
the densities of the injected particles at the fax left and right 
sides, and axe controlled by the voltage source. The shaded 
area represents a more general asymmetric impurity poten- 
tial, which would have a similar asymmetric density profile. 
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2"=* order: 



phase 


exponent z 


Tic ris rric rUs 


(A) 
(B) 


85c 


1 1-1-1 

2 0-2 



3"'* order: 



phase 


exponent z 


Tic ris rric rris Ic h 


(C) 
(D) 


65c + 25s - 2 
24ge - 2 


11 1-1-2 
2 2 -4 



TABLE I: Most relevant exponents z for the rectification cur- 
rent (see Eqs. (1181 and 121ll 'l. Permutations and sign changes 
of n, m, / leading to the same exponents are not shown. 



FIG. 3: Diagram for the leading power-law contributions 
to the rectification current, Ir V^, for voltages V ^ V* . 
The exponents z for the four different phases are given by the 
formulae in the figure (see also Table|TJ. The zone boundaries 
indicate crossover regions. On the left of the dashed line V > 
0, and on the right V — (see Eq. Ull^ 'l. In the hatched 
area, z < 0, and the rectification current shows the behavior 
sketched in Fig. |51 



245c - 2 




2 



